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Abstract 

By using a Generalized Hubbard model for bosons, the energy transfer in a nonlinear 
quantum lattice is studied, with special emphasis on the interplay between local and 
nonlocal nonlinearity. For a strong local nonlinearity, it is shown that the creation 
of v quanta on one site excites a soliton band formed by bound states involving 
v quanta trapped on the same site. The energy is first localized on the excited 
site over a significant timescale and then slowly delocalizes along the lattice. As 
when increasing the nonlocal nonlinearity, a faster dynamics occurs and the energy 
propagates more rapidly along the lattice. Nevertheless, the larger is the number of 
quanta, the slower is the dynamics. However, it is shown that when the nonlocal 
nonlinearity reaches a critical value, the lattice suddenly supports a very fast energy 
propagation whose dynamics is almost independent on the number of quanta. The 
energy is transfered by specific bound states formed by the superimposition of states 
involving v — p quanta trapped on one site and p quanta trapped on the nearest 
neighbour sites, with p = 0, ..,v — 1. These bound states behave as independent 
quanta and they exhibit a dynamics which is insensitive to the nonlinearity and 
controlled by the single quantum hopping constant. 
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1 Introduction 



The concept of energy localization due to nonlinearity in classical lattices has 
been a central topic of intense research over the last three decades. This topic 
can be traced back to the seminal works of Davydov, where the nonlinearity 
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was introduced to explain the vibrational energy flow in a-helices [1]. The main 
idea was that the energy released by the hydrolysis of adenosine triphophate 
(ATP), partially stored in the high-frequency amide-I vibration of a peptide 
group, delocalizes along the helix leading to the formation of vibrational ex- 
citons. Due to their interaction with the phonons of the helix, the excitons 
experience a nonlinear dynamics. Therefore, they propagate according as a 
"Davydov" soliton, a solution of the Nonlinear Schrodinger (NLS) equation 
within the continuum approximation [2,3]. In the mid-1980's, lattice effects 
were introduced through the analysis of the discrete version of NLS. This 
equation supports specific solutions called "lattice solitons" [4] and it has 
revealed the occurrence of a remarkable feature known as the self-trapping 
mechanism [5]. As discovered by Sievers and Takeno [6], the self-trapping is a 
special example of a more general solutions called discrete breathers [7,8,9]. In 
classical anharmonic lattices, a discrete breathers correspond to time-periodic 
and spatially localized solutions which result from the interplay between the 
discreteness and the nonlinearity. These solutions do not require integrability 
for their existence and stability and it has been suggested that they should 
correspond to quite general and robust solutions. Since discrete breathers sus- 
tain a local accumulation of the vibrational energy, which might be pinned in 
the lattice or may travel through it, they are expected to be of fundamental 
importance. 

At present, because the occurrence of classical breathers is a relatively well 
understood phenomena, great attention has been paid to characterize their 
quantum equivalent, for which less detailed results are known [10]. In the 
quantum regime, the Bloch theorem applies due to the translational invari- 
ance of the lattice. Therefore, the corresponding eigenstates cannot localize 
the energy because they must share the symmetry of the translation opera- 
tor which commutes with the lattice Hamiltonian. Nevertheless, the nonlin- 
earity is responsible for the occurrence of specific states called multi-quanta 
bound states [10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27]. A bound 
state corresponds to the trapping of several quanta over only a few neighbour- 
ing sites, with a resulting energy which is less than the energy of quanta lying 
far apart. The distance separating the quanta is small, so that they behave as 
a single particle delocalized along the lattice with a well-defined momentum. 
Since the occurrence of these bound states results from the nonlinearity, they 
can be viewed as the quantum counterpart of breathers or solitons. In low- 
dimensional lattices, two-quanta bound states have been observed in molecular 
adsorbates such as H/Si(lll) [28,29], H/C(lll) [30], CO/NaCl(100) [31] and 
CO/Ru(001) [32,33,34,35,36] using optical probes. Bound states in the system 
H/Ni(lll) were investigated via high resolution electron energy loss spec- 
troscopy [37]. Moreover, a recent experiment, based on femtosecond infrared 
pump-probe spectroscopy, has clearly established the existence of bound states 
in a- helices [38]. 



2 



From a theoretical point of view, most of the previous work was performed 
within the quantum equivalent of the discrete NLS equation. The correspond- 
ing Hamiltonian is essentially a Bose version of the Hubbard model which 
has been used to study a great variety of situations ranging from molecular 
lattices [11] to Bose-Einstein condensates [39,40,41]. Within this model, the 
nonlinearity is local so that it is responsible for a strong interaction between 
quanta located on the same site. However, through the small polaron descrip- 
tion of the vibrational energy flow in proteins, it has been shown recently 
that a nonlocal nonlinearity strongly modifies the nature of the bound states 
[22,23,24,25,42]. In such systems, the polaron results from the dressing of a 
vibrational exciton by a virtual cloud of phonons. In addition to an attractive 
coupling between polarons located onto the same site, a coupling takes place 
when polarons lie on nearest neighbour sites, due to the overlap between their 
virtual cloud of phonons. When two quanta are excited, the competition be- 
tween the local and the nonlocal nonlinearity favours the occurrence of two 
kinds of bound states whose properties have been discussed in details in Ref. 
[22]. 

At this step, the fundamental question arises whether the interplay between 
the local and the nonlocal nonlinearity modifies the dynamics of bound states 
involving several quanta. This is the purpose of the present paper. It will 
be shown, using a generalized Hubbard model, that a critical value of the 
nonlocal nonlinearity favours a resonance responsible for a fast energy transfer. 
Note that, although such an approach is rather general, it will be applied to 
vibrational excitons moving in a nonlinear molecular lattice. 

The paper is organized as follows. In Section 2, the generalized Hubbard model 
is described and the number state method [17], used to define its eigenstates, is 
summarized. Then, the quantum dynamics required to characterize the trans- 
port properties is introduced. The corresponding time dependent Schrodinger 
equation is solved numerically in Section 3 where a detailed analysis of the 
multi-quanta dynamics is performed. Finally, the numerical results are dis- 
cussed and interpreted in Section 4. 



2 Model Hamiltonian and Quantum analysis 

2.1 Hamiltonian and Quantum states 

We consider a one-dimensional lattice formed by N sites whose position is 
defined in terms of the integer index n. Each site n is occupied by a high- 
frequency oscillator described by the standard boson operators and b n . The 
lattice dynamics is governed by a generalized Hubbard model for bosons whose 
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Hamiltonian is written as (using the convention h = 1) 

H = J2 uob% - Abib{b n b n - Bbi +1 b% +1 b n + 6„+i + &l +1 6 n ] (1) 

n 



where o;o is the internal frequency of each oscillator, $ is the hopping con- 
stant between nearest neighbour sites, and A and B represent the local and 
the nonlocal nonlinearity, respectively. Note that for vibrational excitons, the 
nonlinear parameters are positive. 

Since the Hamiltonian H (Eq.(l)) conserves the number of quanta, its eigen- 
states can be determined by using the number state method detailed in Ref. 
[17]. This method is summarized as follows. Because H commutes with the 
operator J2 n Hfin which counts the total number of quanta, the Hilbert space 
E is written as the tensor product E = Eq <g) E\ ® E2 <S> ■■■ <S> E v ..., where 
E v refers to the v -quanta subspace. To generate E v , a useful basis set is 
formed by the local vectors \pi,...,p N ) (J2 n Pn = v), where p n is the num- 
ber of quanta located on the nth site. The dimension d v of E v is equal to the 
number of ways of distributing v indistinguishable quanta among N sites, i.e. 
d v = (N + v- l)\/v\(N- 1)!. 

Within this representation, the Hamiltonian is block-diagonal. However, the 
size of each block can be reduced by taking advantage of the lattice periodicity. 
Indeed, due to the translational invariance, the lattice momentum k is a good 
quantum number so that an "momentum" basis set can be formed. For each 
k value, an momentum vector |<& a (fc)) is obtained by superimposing all the 
vectors \p±, ...,Pn) describing similar quanta distributions related to each other 
by a translation along the lattice. For instance, the momentum vector built 
from the N vectors \v, 0, 0), ... , |0, 0, v ) is defined as 

1 JV-l 

!$!(*;)> = ^]Te te T>,0,...,0> (2) 



where the translation operator T satisfies T\p 1 ,p 2 ...,Pn) = \pn,Pi---,Pn-i)- 

Consequently, in each subspace E v , the Hamiltonian exhibits independent 
blocks associated to each k value and the corresponding Schrodinger equa- 
tion can be solved with a minimum of computational effort. This numerical 
procedure yields the eigenvalues u^\k), which define the dispersion relations 
of the v -quanta subspace, and the associated eigenvectors ^\k). 
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2.2 Quantum dynamics 



The quantum dynamics is governed by the time dependent Schrodinger equa- 
tion written as 

^ = «!•(*)> (3) 



where \^(t)) is the lattice quantum state at time t. By assuming that the state 
is known at the initial time t — 0, its value at time t > is expressed in terms 
of the eigenstates of H as 

!*(*)> = EE^WI^ )) 6 "' 4 '^^?^)) ( 4 ) 



The time dependent quantum state is the central object of the present 

study. Its knowledge allows us to characterize the multi-quanta dynamics 
through the computation of the expectation value of any relevant observable 
O as (0(t)) = (ty(t)\0\ty(t)). Nevertheless, the solution of Eq.(4) requires the 
specification of an initial quantum state. To proceed, we restrict our attention 
to a spatially localized state corresponding to the creation of v quanta on the 
n th site as 

|¥(0)> = -$=|0) (5) 



where |0) denotes the vacuum with zero quanta. This specific choice allows 
us to characterize the ability of the nonlinear quantum lattice to localize the 
energy, at least over a given timescale. Indeed, for a vanishing nonlocal non- 
linearity (B = 0), each subspace E v supports a low energy band called the 
soliton band [16,17]. For vibrational excitons, the nonlinearity A is usually 
greater than the hopping constant $. Consequently, the soliton band describes 
v-quanta bound states in which the v quanta are trapped on the same site and 
behave as a single particle. When v quanta are created on a given state, only 
the soliton band is significantly excited. Because the dispersion of the band 
scales as <& v /A v ~ l when A ^> <£>, the v quanta are localized in the vicinity of 
the excited site over a timescale which increases with both the local nonlinear- 
ity and the number of quanta. This localized behaviour, which is the quantum 
signature of the classical self-trapping, disappears in the long time limit due 
to the non vanishing dispersion of the soliton band. 

In that context, the aim of the present paper is to analyze the way this scenario 
is modified when the nonlocal nonlinearity is turned on. 
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Fig. 1. Energy spectrum 
for N = 91, v = 3 
and A = 3<5 and for (a) 
B = 0, (b) B = 3$ and 
(c) B = m. The left-hand 
side of the figure repre- 
sents the corresponding lo- 
cal density of state (see the 
text). 
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3 Numerical results 



In this section, the numerical diagonalization of the Hamiltonian H is per- 
formed and the time dependent Schrodinger equation is solved. To realize the 
simulation, the hopping constant $ is used as a reduced unit and the local 
nonlinearity is fixed to A = 3$. The nonlocal nonlinearity is taken as a free 
positive parameter. 

When v — 3, the energy spectrum of the Hamiltonian H is displayed in Fig. 
1 for B = (Fig. la), B = 3$ (Fig. lb) and 5 = 6$ (Fig. lc). In each 
figure, the spectrum is centred onto the frequency 3ujq and it corresponds to 
the dispersion curves drawn in half of the first Brillouin zone of the lattice, 
i.e. < k < ix. The left-hand side of Fig. 1 represents the corresponding local 
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density of state (LDOS). It describes the weight of the initial state |^(0)) in 
the eigenstates the energy of which ranges between ui and uj + du as 

PoH = 0*i 3) WI*(°)}| 2 % - (6) 

kX 

When B = (Fig. la), the energy spectrum exhibits two energy continua 
and three isolated bands. The high frequency continuum characterizes free 
states describing three independent quanta. The low frequency continuum 
supports states in which two quanta are trapped on the same site whereas 
the third quantum propagates independently. By contrast, the isolated bands 
refers to bound states. The low frequency band is the well known soliton band. 
It lies below the continua over the entire Brillouin zone and it characterizes 
three quanta trapped on the same site and delocalized along the lattice. The 
two other bands, which lie respectively below and above the low frequency 
continuum, describe bound states in which two quanta are trapped on a given 
site, whereas the third quantum is trapped onto the corresponding nearest 
neighbour sites. As shown in the left-hand side of Fig. la, the LDOS is strongly 
peaked in the frequency range of the soliton band, which is significantly excited 
when three quanta are initially created on the same site. 

When 5 = 3$ (Fig. lb), the energy spectrum supports three continua and 
four isolated bands. Indeed, a third continuum occurs between the two previ- 
ous energy continua. It describes states for which two quanta are trapped on 
two nearest neighbour sites whereas the third quantum propagates indepen- 
dently. The low frequency isolated band is the soliton band, whose width has 
been increased when compared with the previous situation. Above the soli- 
ton band, the next two isolated bands correspond to the previously observed 
bound states in which two quanta are trapped on one site whereas the third 
quantum is localized on nearest neighbour sites. Note that these bands have 
been strongly redshifted. Finally, a fourth bound state band occurs just below 
the continua. The analysis of the corresponding eigenstates reveals that this 
band supports bound states involving three quanta trapped on three nearest 
neighbour sites. In addition to the occurrence of this latter band, the behaviour 
of the LDOS reveals that the non vanishing B value is responsible for a weak 
hybridization between the three lowest bound state bands. In other words, 
although the soliton band still mainly refers to three quanta on the same site, 
the two other bound states involve such a configuration but to a lesser extent. 

Finally, when B = 6$ (Fig. lc), the energy spectrum still supports three con- 
tinua and four isolated bands. The shape of the continua has been slightly 
modified and the high frequency bound state band, which still mainly refers 
to three quanta trapped onto three neighbouring sites, has been slightly red- 
shifted. However, strong modifications affect the behaviour of the three low 
energy bound state bands. Indeed, the gaps between these bands are now very 
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Fig. 2. Space and time evolution of the number of quanta for N = 91 and v = 3 and 
for (a) B = 0, (b) B = 2$, (c) B = 4$, (d) B = 6$, (e) B = 8$ and (f) B = 10$ 

small, which indicates the occurrence of a strong hybridization. The LDOS 
exhibits a significant value over the frequency range of the three low energy 
bands. It shows four peaks for specific frequencies, corresponding either to the 
bottom or to the top of the three low frequency bands (see Fig. lc). In other 
words, these bands support bound states which are superpositions of states 
involving three quanta trapped on one site and two quanta on one site trapped 
with a third quantum localized on neighbouring sites. 

In Fig. 2, the time evolution of the expectation value of the population op- 
erator, P n (t) = (^ r (t)|6 n 6J l |^ r (i)), is displayed for v — 3 and for different B 
values. The size of the lattice if fixed to N — 91 sites. When B = (Fig. 
2a), the population of the initial site P no (t) decreases very slowly with time. It 
reaches 50 % of its initial value when t is about 13. 5$ -1 , which indicates that 
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Fig. 3. Space and time evo- 
lution of the number of 
quanta for N = 29 and 
B = 2A and for (a) v = 2, 
(b) v = 3 and (c) v = 4. 



the energy stays localized on the excited site no over a significant timescale. 
Nevertheless, energy propagation occurs and the population delocalizes along 
the lattice. However, this propagation corresponds to a rather slow process, 
since 95 % of the popuplation is concentrated in the 6 sites surrounding no 
at time t = SO^" 1 . When B = 2$ (Fig. 2b), the same features occur. Never- 
theless, a faster dynamics takes place since P no {t) reaches 50 % of its initial 
value for t = 6.8$ _1 . The energy propagates slowly and 95 % of the popu- 
lation is concentrated in the 10 sites surrounding n for t = 30$ _1 . When 
B = 4$ (Fig. 2c), the population of the excited site still decays but it now 
supports high frequency oscillations. Its decay is enhanced, and it reaches 50 
% of its initial value when t = 2.1$~ 1 . The population propagates more easily 
along the lattice and its maximum value is observed on the sites no ± 12 for 
t = 30$ _1 . In that case, 95 % of the population is contained in the 30 sites 
surrounding no- 
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When 5 = 6$ (Fig. 2d), a fully different behaviour is observed. Indeed, P no (t) 
shows a very fast decay since it reaches 50 % of its initial value for t = 0.6$ _1 . 
This decay occurs over a timescale of about 20 times shorter than the time 
associated to the decay of the population for 5 = 0. In addition, Fig. 2d 
clearly shows that the lattice supports a rather fast energy propagation. Two 
population wave packets are emitted on each side of the initial site with a 
velocity of about $. In other words, the maximum value of the population is 
observed on the sites n ± 35 for t = 30$ _1 . 

Finally, for greater 5 values, a rather slow dynamics recurs. When 5 = 8$ 
(Fig. 2e) and B = 10$ (Fig. 2f), the population of the excited site slowly 
decreases. It supports a small amplitude high-frequency modulation, and it 
reaches 50 % of its initial value for t = 4.2$ _1 and t = 6.3$~ x when 5 = 8$ 
and B = 10$, respectively. The population propagates along the lattice with 
a rather small velocity. Therefore, when t = 30$ _1 , 95 % of the energy is 
concentrated on the 18 sites surrounding the excited site when B = 8$. The 
size of this region decreases to 14 sites when B = 10$. 

Fig. 2 has clearly shown the occurrence of a specific dynamics when the relation 
B = 2A is verified. As illustrated in Fig. 3, such a behaviour does not depend 
significantly on the value of the number of quanta. Indeed, the figure shows 
the time evolution of the population for v — 2 (Fig. 3a), v — 3 (Fig. 3b) and 
v = 4 (Fig. 3c) and for B = 2A. Note that the lattice size is fixed to N = 29 
to avoid an overlong computational time, especially when v = 4. Whatever 
the number of quanta, the population of the excited site decays over a short 
timescale. Then, two population wave packets are emitted on each side of the 
excited site with a velocity which slightly decreases with the number of quanta. 
When v — 2, the wave packets reach the edges of the lattice for t = 10.8$ _1 
whereas longer delays t = 12.6$ _1 and t = 13.7$~ 1 are required when v — 3 
and v = 4, respectively. 

The time evolution of the survival probability is displayed in Fig. 4 for v — 2 
(Fig. 4a), v — 3 (Fig. 4b) and v — 4 (Fig. 4c). The number of sites is fixed 
to iV = 29 and three different B values have been used, i.e. 5 = (full line), 
5 = 3$ (dashed line) and 5 = 6$ (dotted line). The survival probability 
So(t) is defined as the probability to observe the lattice in the state |^(0)) 
at time t, i.e. S (t) = \(^(0)\^(t))\ 2 . It is expressed in terms of the Fourier 
transform of the LDOS (Eq.(6)) as 



When 5 = 0, the survival probability slowly decreases with time, which indi- 
cates that the energy is localized on the excited site over a significant timescale 




(7) 



— oo 
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which increases with v. For instance, the first zero of So(t) is reached for 
t = 4$ _1 and t = 30$ _1 when v — 2 and 3, respectively. When v — 4, So(t) 
is almost constant over the timescale displayed in the figure. Note that So(t) 
supports a small amplitude high-frequency modulation which is clearly seen 
for v = 3 and v — 4. As when increasing B, S (t) exhibits a faster decay. 
When B = 3$, the first zero of the survival probability occurs for t = 2.4$ _1 
and t = 9.4$ _1 when v — 2 and 3, respectively. In a similar way, although the 
first zero cannot be observed for v — 4, S'o(t) decreases with time (Fig. 4c). In 
a marked contrast with the case B = 0, the amplitude of the modulation in- 
creases with B. Finally, when B = 6$, S (t) suddenly shows a very fast decay 
whatever the number of quanta. This decay occurs over a similar timescale 
and a first minimum is reached for t = 0.9$ _1 , t = 0.75$~ 1 and t = 0.7$ _1 
when v — 2, 3 and 4, respectively. In addition, due to the small lattice size, 
a revival is observed for t = 20$ _1 , t = 25$~ 1 and t = 27$~ 1 when v = 2, 3 
and 4, respectively. These latter times represent the delay for the excitation to 
cover the whole lattice, so that the corresponding velocities are about 1.44$, 
1.16$ and 1.06$. Therefore, in a perfect agreement with the results observed 
in Fig. 3, this feature indicates that when B = 2A a fast energy propagation 
occurs. The initial excitation covers the lattice over a timescale almost inde- 
pendent on the number of quanta and with a velocity typically of the order of 
the velocity of a single quantum insensitive to any nonlinearity. 

To determine the relevant number states which contribute significantly to the 
quantum dynamics, we introduce W a (t) as the sum over k of the weight of 
each intermediate basis vector |$ Q (A;)) as 

w«(0 = EK $ «(*)l*(*)>r (8) 

k 

Since a number state refers to a given distribution of the quanta, the asso- 
ciated momentum vector characterizes all the local states connected to that 
distribution but related to each other through a translation along the lattice. 
Therefore, the sum over k allows us to characterize the participation of such a 
distribution in the quantum dynamics. For instance, from the N local vectors 
\v, 0, 0), ... , |0, 0, v), we have defined in Eq.(2) the vector |$i(fc)) which 
basically describes v quanta trapped on the same site. The corresponding 
weight Wi(t) allows us to quantify the participation of such a configuration. 

For v — 3, Fig. 5 displays the relevant number states involved in the quantum 
dynamics for B = (Fig. 5a), B = 3$ (Fig. 5b) and 5 = 6$ (Fig. 5c). 
When B = 0, the states involving 3 quanta trapped on the same site control 
more than 90 % of the quantum dynamics. However, as when increasing B, 
the participation of other states takes place. When 5 = 3$ (Fig. 5b), the 
states describing 3 quanta trapped on the same site remain the most relevant 
since they carry about 75 % of the full dynamics. However, about 22 % of the 
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Fig. 4. Survival probabil- 
ity for N = 29 and for 
(a) v = 2, (b) v = 3 and 
(c) -u = 4. Three -B values 
have been used, i.e. B = 
(full line), 5 = 3$ (dashed 
line) and B = 6® (dotted 
line) . 



dynamics is controlled by states involving two quanta on one site trapped to 
a third quantum located on the nearest neighbour site. A factor 2 is included, 
due to the participation of the two equivalent states 1 2100. . . ) and 1 1200...) . 
Finally, when 5 = 6$ (Fig. 5c), the dynamics are mainly controlled by the 
states referring to two quanta on one site trapped to a third quantum. These 
states, including the factor two, carry about 56 % of the full dynamics whereas 
35 % originates in states involving 3 quanta trapped on one site. Note that 
a third type of states connected to three quanta trapped onto three nearest 
neighbour sites participate to a lesser extent, since they share only 8 % of the 
dynamics. 

Finally, the different observed features are illustrated in Fig. 6 for v — 5. The 
lattice size is fixed to iV = 19 to avoid an overlong computational time. Nev- 
ertheless, Fig. 6a shows the occurrence of the derealization of the population 
when B = 2A. As when v = 2, 3, and 4, a fast decay of the population of 
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Fig. 5. Relevant num- 
ber states involved in the 
quantum dynamics for (a) 
B = 0, (b) B = 3$ and (c) 
B = 6$. 



the excited site takes place. It is followed by the emission of two wave packets 
which propagate according to velocity of about <£>. In Fig. 6b, the survival 
probability indicates a strong energy localization when B = and B = 3$, 
whereas it decreases rapidly over a timescale of about 0.5$ _1 when B = 2A. 
Finally, Fig. 6c shows that the main part of the quantum dynamics is described 
by the states involving the configurations 1 500. . .) (about 26 %), |410...) and 
|140...) (about 40 %), |320...) and |230...) (about 22 %) and |131...) (about 8 
%). 
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Fig. 6. Lattice dynamics 
for v = 5 and N = 19. 
(a) Space and time evo- 
lution of the number of 
quanta, (b) survival prob- 
ability and (c) relevant 
number states. 
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4 Interpretation and discussion 



In the previous section, the numerical results have revealed a strong depen- 
dence of the lattice dynamics on the nonlocal nonlinearity B. Indeed, for a 
vanishing B value, the creation of v quanta on one site mainly excites the 
soliton band of the energy spectrum. This excitation gives rise to the occur- 
rence of the quantum equivalent of the classical self-trapping phenomena. The 
energy is localized on the excited site over a significant timescale, which in- 
creases with both the local nonlinearity and the number of quanta. Then, a 
rather slow energy flow takes place due to the finite value of the soliton band- 
width. As when increasing the nonlocal nonlinearity, a faster dynamics occurs. 
The survival probability, which characterizes the memory of the initial state 



14 



in the time dependent quantum state, decreases more rapidly and supports 
a significant high frequency modulation. Similarly, the population propagates 
more rapidly along the lattice. In that regime, the speed of the dynamics still 
depends on the number of quanta and the larger is the number of quanta, 
the slower is the dynamics. However, when the nonlocal nonlinearity reaches 
the critical value B = 2A, a very different behaviour takes place. Indeed, the 
lattice supports a very fast energy propagation whose dynamics is almost in- 
dependent of the number of quanta. The survival probability decreases over a 
timescale of about and two population wave packets emitted on each side 
of the excited site propagate with a velocity typically of about $. Our numer- 
ical results have clearly established that this fast energy transfer is mediated 
by bound states. However, the timescale of the dynamics is of the same or- 
der of magnitude as the timescale governing the free state dynamics. In other 
words, the bound states involved in the dynamics when B = 2A behave like 
independent quanta insensitive to the nonlinearity. Finally, when B exceeds 
its critical value, the quantum self-trapping regime recurs. 

As shown in Section 3, the main part of the dynamics is controlled by specific 
number states. Indeed, we have verified numerically for v = 2,..., 7 that the 
most relevant states involve number states of the form \v — p,p, 0, 0), with 
p — 0, .., v — 1. Note that the states |1, v — 2, 1, 0), |1, v — 3, 2, 0) and 
\2,v — 3, 1,...,0) participate to a lesser extend in the dynamics. Therefore, 
the relevant states, which describe v — p quanta on one site trapped with 
p quanta on a nearest neighbour site, are characterized by the self-energies 
e p = vu -v(v - l)A + p(v - p)(2A - B) (see Eq.(l)). As long as \2A - B\ is 
strong enough, the different relevant configurations are weakly coupled to each 
other. Therefore only the number states |i>, 0, 0, 0) participate significantly 
to the dynamics which follows the initial creation of v quanta on one site. 
However, as illustrated in Fig. 5, a strong hybridization between the different 
configurations occurs when \2A — B\ tends to zero. Finally, when B = 2 A, a 
resonance takes place since all the relevant states have the same energy. 

To understand more clearly the way this resonance modifies the energy trans- 
fer in the nonlinear lattice, we can take advantage of the fact that only relevant 
number states participate in the full dynamics. This feature allows us to es- 
tablish a simplified model, which is able to account for the dynamics in the 
case A>$. This model is introduced in the next section. 



4-1 Equivalent lattice model 

As mentioned previously, a good description of the dynamics in the subspace 
E v is obtained by restricting the number state basis to the set of the relevant 
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Fig. 7. Equivalent lattice 
model. 



\n - 1,0) 



\n,0) \n,l) \n,2) |n,3> \n + 1,0) 



(b) 



|n-l> 



|n) 



|n + l> 



vectors. This set is formed by the N x v vectors defined as 



n,p) = \pi = 0,p 2 = 0,...,p n = v-p,p n+1 =p,...0) 



(9) 



where \n,p) characterizes v — p quanta on the nth site and p quanta of the 
site n + 1. The representation of the Hamiltonian H in the restricted basis 
is equivalent to a tight-binding model on the lattice displayed in Fig. 7a. 
To simplify the notation, the restriction of the full Hamiltonian will still be 
denoted H. Within this model, each site supports the state \n,p) whose energy 
is e p = p(v — p)e, where e = 2A — B. Note that eo = has been used as 
the origin of the energy. Nearest neighbour sites are coupled to each other 
through generalized hopping constants. The hopping constant between \n,p) 



and \n,p + 1) is equal to $ p = y (p + l)(f — with p = 0, .., v — 2. By 
symmetry, the hopping constant between \n, v — 1) and \n + 1, 0) is <3> = y 7 ^- 

Consequently, the dynamics of the quantum nonlinear lattice is formally equiv- 
alent to the dynamics of a fictitious particle moving on the lattice shown in 
Fig. 7a. At time t = 0, the particle is created in the state |n , 0) and its time 
evolution is governed by the corresponding Schrodinger equation. Unfortu- 
nately, although the model has been greatly simplified, this equation cannot 
be solved analytically. However, the knowledge of some general properties is 
sufficient to clarify our understanding of the physics involved in the energy 
transfer. In addition, the model allows us to determine the LDOS which char- 
acterizes the weight of a localized state |n, 0) in the eigenstates of the lattice. 
It represents a measure of the response of the lattice to the initial excitation 
and gives the survival probability simply by performing a Fourier transform. 

4-2 General properties of the lattice eigenstates 

The lattice displayed in Fig. 7a exhibits a translational invariance and its 
unit cell is defined in terms of two subsystems. The first subsystem is formed 
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by the single state \n, 0), whereas the second subsystem involves the v — 1 
coupled states \n,p), with p — 1, .., v — 1. A moment's reflection will convince 
the reader that this latter subsystem represents a modified nonlinear quantum 
dimer. More precisely, it is the representation of a dimer in the number state 
basis in which the two states involving v quanta on the first site and v quanta 
on the second site have been suppressed. As for the dimer, this subsystem is 
invariant under the reflection symmetry which transforms the state \n,p) into 
the state \n, v — p). 

Due to the translational invariance, the lattice wave vector A; is a good quan- 
tum number. Therefore, for each k value, the quantum states are described 
by a v x v matrix H(k). Instead of representing this matrix in the local basis 
\n,p), it is more convenient to introduce a new basis set formed by the sin- 
gle state | ra, 0) and by the v — 1 eigenstates of the modified dimer. Indeed, 
although the quantum states of a modified dimer are only known for small v 
values, they represent a key ingredient needed to understand the physics of 
the lattice. Therefore, let \a) denotes the ath eigenstate of a modified dimer 
whose energy is oj a (see Fig. 7b). Because of the reflection symmetry, each 
eigenstate is either symmetric or antisymmetric. As a result, the wave func- 
tion tp ap = (n,p\a) satisfies ip ap = ±ip av _ p , depending on whether the state is 
symmetric (+) or antisymmetric (— ). For even v values, there are N$ — v/2 
symmetric and Na = (v — 2)/2 antisymmetric eigenstates. By contrast, for 
odd v values, N s = N A = (v - l)/2. 

In that context, the matrix Hamiltonian H{k) is expressed as 



H(k) 



' x a (k) x a f(k) ... N 
x* a {k) uj a .. 



x* a ,(k) 



V 



Uryt 



(10) 



/ 



where x a (k) = y/v(p{ip a i + ip a v-i exp(— ik)) is the coupling between the state 
\n, 0) and \a) which results from the interaction between \n, 0) with both the 
state | n, 1) of the same unit cell and the state \n — l,v — 1) of the previous 
unit cell (see Fig. 7). 

iFrom Eq.(10), the derealization of the fictitious particle originates in the 
coupling between two neighbouring states \n, 0) and \n+ 1,0), through their 
interaction with the eigenstates of the nth modified dimer. As a consequence, 
each lattice eigenstate is a superimposition involving the localized states and 
the eigenstates of a modified dimer. However, such a hybridization depends 
both on the energy difference cu a , and on the strength of the coupling x a (k), 
which is drastically sensitive to the symmetry of the dimer eigenstates. 
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Eq.(lO) allows us to determine the exact expression of the LDOS. Indeed, 
according to the standard definition used in condensed matter physics, the 
LDOS is expressed in terms of the Green's operator G(k,cu) = {00 — H(k))^ 1 
as 



Po(u) = 




Since the LDOS only depends on the restriction of the Green's operator to 
\no, 0), it can be obtained by applying standard projection methods. Therefore, 
by taking advantage of the symmetry of the eigenstates \a), the restricted 
Green's operator is written as 

(n o ,0\G(k,cu)\n o ,0) = 1 (12) 

to — A(u) — 2J(u) cos(fc) 



where 



A( W ) = 2# 2 ^ ilil 



CO — LO r 



on — 1 /, \ 

J(w) = v$ 2^— — (13) 

a LO — LO a 

By inserting Eq.(12) into Eq.(ll), the LDOS is finally written as 

PoH = - 1 (14) 

^ ^{2J{uo) + A(w) - oo){2J{u) - A(w) + to)) 



Eq.(14) allows us to establish two general properties. First, the zeros of the 
LDOS are given by the poles of the parameters A (a;) and J(oo). In that context, 
Eq.(13) clearly shows that the LDOS supports v — 1 zeros, which corresponds 
to the eigenenergies of the modified dimer. In the vicinity of a zero, the LDOS 
scales as \J\u — u a \. Then, the LDOS exhibits two kinds of poles which are 
given by the solutions of the two equations 2J(uo) ± (A(cj) — uo) = 0. From 
Eq.(13), these two equations are rewritten as 



S 12 
A / 2 

Av &yJ^= u (15) 

where the symbol J2 S,A denotes a sum over either the symmetric or the anti- 
symmetric eigenstates of a modified dimer. In Eq.(15), the first equation gives 
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Ns + 1 solutions whereas the second equation yields Na + 1 solutions. There- 
fore, the LDOS supports Ng + Na + 2 = v + 1 poles. In the vicinity of a pole 
Upo, the LDOS thus diverges according to the power law \uj — uj po \~ l l 2 . 

These results reveal the occurrence of zeros and poles in the LDOS, as observed 
in Fig. 1. In addition, they point out the important role of the symmetry of 
the eigenstates of the dimer, especially for specific values of the lattice wave 
vector k. Indeed, when k = 0, the coupling x a (k = 0) between \n, 0) and each 
antisymmetric state \a) vanishes. As a consequence, H(k = 0) (Eq.(10)) has 
Na eigenvalues equal to the eigenenergies of the antisymmetric states of the 
modified dimer. They describe lattice eigenstates which do not involve |n, 0), 
so they corresponds to zeros of the LDOS. By contrast, the remaining coupling 
between \n, 0) and the symmetric states \a) produces N$ + l eigenstates which 
induce divergences in the LDOS. When k = n, the same features occur, but 
by inverting the role played by the symmetric and the antisymmetric states. 
The matrix H(k = n) exhibits N s eigenvalues equal to the eigenenergies of the 
symmetric states \a) and which correspond to zeros of the LDOS. A strong 
coupling remains between \n, 0) and the Na antisymmetric states, leading to 
Na + 1 eigenvalues which give rise to divergences in the LDOS. 

Finally, the coupling between \n, 0) and the v — 1 dimer eigenstates \a) is at 
the origin of the creation of the lattice eigenstates. It is responsible for the 
occurrence of specific signatures in the LDOS and control the lattice dynam- 
ics. However, depending on the strength of the coupling, different dynamical 
behaviours can be observed ranging from quantum self-trapping to fast energy 
transfer. In the following of the text, these features are first illustrated for the 
two simples situations corresponding to v — 2 and v = 3. Then, a general 
discussion is given to interpret the specific behaviour observed when e = 0. 



4-3 Application to v = 2 and v = 3 



The case v — 2 is rather simple, because a modified dimer involves a single 
state |n, 1) whose energy is e\ = e. Therefore, H(k) is a 2 x 2 matrix which 
describes the hybridization between \n, 0) and \n, 1). It can be solved exactly 
so that the lattice supports two bands whose dispersion relations are expressed 
as 




(16) 
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The corresponding LDOS, which shows a single zero and three poles, is defined 

as 



1 



-(w-ei) 



(17) 



7T 



When e ^> $ (B <C 2A), the two bands are only weakly coupled. The low 
frequency band, whose energy is about e = 0, refers to the derealization of 
the fictitious particle over the different states |n, 0). By contrast, the second 
band accounts for the derealization of the fictitious particle over the states 
| n, 1). In other words, the low frequency band is the soliton band describing 
two quanta trapped on the same site and delocalized along the lattice. The 
second band refers to the derealization of bound states involving two quanta 
trapped onto two nearest neighbour sites. 

Since only the soliton band is significantly excited, it controls the dynamics. To 
understand this feature, let us evaluate the corresponding survival probability 
(Eq.(7)). Although the Fourier transform of Eq.(17) is not analytic in a general 
way, it can be evaluated when the two bands lie far from each other. In that 
case, the integral over the frequency range of a given band can be performed by 
neglecting the frequency dependence of the other band. After straightforward 
calculations, the survival probability is thus expressed as 



where Jo is the Bessel function of the first kind. The survival probability slowly 
decreases with time and it reaches its first zero for t = 2.4e/4$ 2 . Note that 
this value is extracted from the first zero t 2.40 of Jo(t). When A = 3$ 
and B = 0, this value is equal to 3.6$ _1 , in a perfect agreement with the 
numerical results about 4<3> -1 (see Fig. 4a). Eq.(18) shows that the survival 
probability supports a high frequency modulation whose amplitude decreases 
with e. Such a modulation, whose frequency is equal to the energy difference 
between \n, 0) and \n, 1), characterizes the small participation of the second 
band in the dynamics. 

As when e decreases, i.e. when B tends to 2A, the width of the low frequency 
band increases as a result of the redshift of its low frequency edge. By contrast, 
the second band exhibits a strong redshift accompanied by an increase in 
width. In addition, its contribution to the LDOS increases, which indicates 
that the coupling between the two kinds of bound states is enhanced. Such a 
behaviour occurs until the high frequency edge of the soliton band meets the 
low frequency edge of the second band. This feature takes place when e = 0, 
i.e. when B = 2A. In that case, the LDOS displays a single band, only, with 



4$ 2 t 8$^ 
So(*)« Jo( )[1-— (l-cos(et))] 



(18) 
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two divergences whose frequencies are ±2y / 2 ( 3>. The corresponding bandwidth 
is thus a maximum and it is equal to 4^$. At the resonance, this band 
clearly suggest strong hybridization between \n, 0) and \n, 1). In that case, 
the Fourier transform of Eq.(17) can be determined exactly and it yields the 
following survival probability 

S (t) = Jg(2>/2$f) (19) 



The resonance induces a decay of the survival probability over a very short 
timescale. It reaches its first zero for t = 0.85<I> _1 , in perfect agreement with the 
results shown in Fig. 4a. This feature, mainly due to the strong hybridization 
between the two bands, is responsible for a fast energy transfer along the 
whole lattice. In other words, the transport of energy results from a series of 
transitions between states involving successively two quanta on the same site 
and two quanta on nearest neighbour sites. 

When v — 3, a modified dimer involves the two states \n, 1) and \n, 2), which 
have the same energy 2e. The first state describes two quanta on site n and 
one quantum on site n + 1, whereas the second state refers to one quantum 
on site n and two quanta on site n + 1. Due to the coupling $i = 2$ be- 
tween these states, a strong hybridization takes place so that the eigenstates 
of the modified dimer correspond to a symmetric and to an antisymmetric 
superimposition as 

|a±) = ^(|n,l)±|n,2)) (20) 



The corresponding eigenenergies are u± = 2e ± 2$. Therefore, H(k) is a 3 x 3 
matrix which describes the hybridization between \n, 0) and \a±). In a general 
way, the lattice supports three bands and the corresponding LDOS is written 

as 



-(u - u + )(u - U-) 
(u - E^u - E 2 )(u - E 3 )(u - E 4 ) 1 ' 

where the four poles are defined as 

E! = e-$- Ve 2 - 2e$ + 7$ 2 
E 2 = e + $ - Ve 2 + 2e$ + 7$ 2 
E 3 = e - $ + Ve 2 - 2e$ + 7$ 2 

E 4 = e + $ + Ve 2 + 2e$ + 7$ 2 (22) 



7T\ 
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For strong e values (see Fig. 8a), the three bands are well separated, in a perfect 
agreement with the results displayed in Fig. la and lb. As for v — 2, the low 
frequency band mainly describes the derealization of the fictitious particle 
over the different states |n, 0). It thus corresponds to the soliton band, the 
LDOS of which exhibits two divergences whose frequencies are E\ and Ei- The 
resulting bandwidth is thus about 6$ 3 /e 2 . By contrast, the two other bands 
account for the derealization of the fictitious particle over the eigenstates \a±) 
belonging to different unit cells. Consequently, the band just above the soliton 
band refers to bound states formed by the antisymmetric superimposition of 
the states \n, 1) and \n, 2), whereas the high frequency band characterizes 
bound states involving their symmetric superimposition. As shown in Fig. 8a, 
the LDOS connected to these two bands shows a single divergence and a single 
zero. 



Since only the soliton band is significantly excited when three quanta are 
created on one site, it controls the main part of the dynamics. The survival 
probability is thus approximately expressed as 



3$ 2 X T9/ 3$ 3 t 



s (t) « (i - JlC—r) 



3$ 2 T ,3$ 3 t, r T ,3$ 2 t, . . T ,3$ 2 t, . , n 
+ _-J (— — ) J (— — ) cos(2et) - J x {—-) sm 2et cos(2$t) 23 
e z e z 2e 2e 
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The survival probability exhibits a rather slow dynamics governed by the first 
term in the right-hand side of Eq.(23). It decreases over a timescale typically 
of about t = 2.40e 2 /3<3> 3 . When A = 3$ and B = 0, this time is equal to 
28. 8$ -1 , in a perfect agreement with the numerical results 30$ _1 (see Fig. 
4b). Note that the survival probability supports a high frequency modulation 
whose amplitude decreases with e. Such a modulation, whose frequency is 
typically about tv±, accounts for the small hybridization between the soliton 
band and the two other bands. 

As for the case when e decreases, i.e. when B tends to 2A, the three bands are 
redshifted and their bandwidth increases. In addition, the contribution to the 
LDOS of the two high frequency bands is enhanced due to their hybridiza- 
tion with the soliton band. Nevertheless, the three bands do not hybridize 
in a similar way. Indeed, the soliton band hybridizes first with the antisym- 
metric band. This mechanism dominates when e = $/4, for which both the 
soliton band and the antisymmetric band condense into a single band. Then, 
for smaller e values, the hybridization with the symmetric band takes place. 
Finally, when e = 0, the LDOS exhibits three bands again (Fig. 8b). However, 
all these bands mix the three different kinds of bound states so that the three 
bands control the dynamics. As a result, the typical time governing the decay 
of the survival probability is about t = 2 x 2.4/Aa>, where Acu = 2(1 + \/7)$ is 
the difference between the high frequency and the low frequency divergences 
of the LDOS. When A = 3$ and B = 6$, this time is equal to CGG^ 1 
in perfect agreement with the numerical results 0.75$ _1 (see Fig. 4b). As a 
consequence, a fast energy transfer takes place, with a velocity of about the 
hopping constant $, which characterizes the motion of single quantum insen- 
sitive to the nonlinearity. However, as for v — 2, the transport of energy is 
mediated by states formed by the superimposition of the three kinds of bound 
states. 



4-4 Interpretation of the resonance e = 

The previous examples clearly illustrate the fact that the lattice dynamics 
exhibits two distinct behaviours, depending on the value of the parameter 
e = 2A- B. 

When e ^> <£> (B <C 2A), the energy transfer is mediated by bound states 
formed by v quanta trapped on the same site and behaving as a single particle. 
This particle delocalizes along the lattice due to the coupling between nearest 
neighbour states \n, 0) through their interaction with the modified dimers. 
Because the eigenenergies of the dimer lie far from the energy of the localized 
state | n, 0), this coupling is very weak. Consequently, the v quanta bound 
states belong to the low energy soliton band, which is characterized by a 
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very small bandwidth of about Au; ~ 4v<&" /(v — \)\e v ~ l [16,17]. Therefore, the 
creation of v quanta on one site mainly excites the soliton band, leading to the 
occurrence of the quantum equivalent of the classical self-trapping. The energy 
is localized on the excited site over a significant timescale which increases with 
both e and v. It finally propagates slowly at a very small velocity, proportional 
to the bandwidth. 

By contrast, when e = 0, i.e. when B = 2A, a very different behaviour takes 
place, mainly due to the occurrence of a resonance, since all the states \n,p) 
have the same self-energy. This resonance allows for a complete energy der- 
ealization at a rather large velocity, typically of about the hopping constant 
$. In a marked contrast with the previous situation, the resonance induces a 
strong interaction between the localized states and the modified dimers. Nev- 
ertheless, as detailed below, the origin of this interaction strongly depends on 
whether v is even or odd, due to the v dependence of the energy spectrum of 
a modified dimer. 

For even v values, a modified dimer supports a single eigenstate, denoted as 
\cxo), whose energy vanishes. The v — 2 remaining eigenstates are grouped into 
pairs formed by two eigenstates having the same symmetry but with opposite 
energy. As a result, the coupling between \n, 0) and all these v — 2 eigenstates 
appears extremely weak. Indeed, by using a standard perturbation theory, it 
is easy to show that the second order correction of the energy of \n, 0) due 
to these couplings is expressed as AE = — J2 a ^a \ x a(k)\ 2 /uj a . This correction 
vanishes exactly because the sum involves pairs of states which experience the 
same interaction with \n, 0) but which have an opposite energy. In other words, 
a pair of eigenstates yields two pathways for the transition between \n, 0) 
and \n + 1,0) which interfere destructively so that the resulting probability 
amplitude vanishes. 

Consequently, the relevant part of H(k) reduces to a 2 x 2 matrix which de- 
scribes the hybridization between \n, 0) and the resonant eigenstate \ao). From 
Eq.(10), it can be solved exactly so that the lattice supports two energy bands 
whose dispersion relations are expressed as iy^^lVw + i'aov-i exp(— ik)\. 
Depending on whether \olq) is symmetric or antisymmetric, these two bands 
vanish for k = it or k = 0. Therefore, the gap between these two bands 
vanishes, so that the resulting bandwidth is equal to Au = Ay/v^\tp ao i\. By 
performing the numerical diagonalization of the dimer Hamiltonian for v rang- 
ing from 2 to 200, it is straightforward to show that the wave function behaves 
as \ip ao i\ ~ \J^/ V - The bandwidth is thus of about a few times the hopping 
constant as Auj ps A\/2<&. Note that the diagonalization of the full H(k) ma- 
trix Eq.(10) shows that the bandwidth slightly decreases with the number of 
quanta. For instance it varies from 5.66$ for v — 2 to 4.50$ for v = 50. This 
rather large bandwidth is responsible for the decay of the survival probability 
over a short timescale, which indicates the occurrence of a fast energy transfer 
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at a velocity typically of about the velocity of independent quanta insensitive 
to the nonlinearity. Nevertheless, this derealization originates in the propa- 
gation of bound states, which result from the hybridization between -y-quanta 
trapped on the same site and f-quanta distributed according to the reso- 
nant state \oto). Our numerical analysis reveals that this distribution involves 
the superimposition of the number states \v — 1,1,0,0...), \v — 3,3,0,0...), 
|w-5,5,0,0...) ...\l,v- 1,0,0...). 



For odd v values, the modified dimer has v — 1 eigenstates which are grouped 
into pairs containing two states having an opposite energy and a different sym- 
metry. Therefore, since no resonance occurs with the modified dimer, \n, 0) in- 
teracts with all the eigenstates \a). Nevertheless, among all these interactions, 
only the coupling between \n, 0) and the two eigenstates whose energy is close 
to zero is significant. These two states, denoted by \as) (symmetric) and \a A ) 
(antisymmetric) respectively, belong to the same group and have an opposite 
energy us = —uj a . Consequently, the relevant part of the H{k) reduces to a 
3x3 matrix written as 



H(k) 



x s (k) x A (k) 
x* s (k) lu s 
x* A (k) -u s 



\ 



(24) 



where x s (k) = \/v^tpsi(^ + exp(— ik)) and x A (k) = v^^si^l — exp(— ik)). 
As discussed in Section 4.2, Eq.(24) can be diagonalized for specific values of 
the lattice wave vector k by taking advantage of the dimer eigenstate sym- 
metry. Indeed, when k = 0, the coupling between \n, 0) and \a A ) vanishes. 
The energy uj a is thus an eigenvalue of the matrix Hamiltonian Eq. (24) which 
corresponds to a zero of the LDOS. However, the coupling between \n, 0) and 

\as) remains and it produces two eigenvalues ws/2 ± ^ (us/2) 2 + 4v^ 2 \ipsi\ 2 
which correspond to two poles of the LDOS. In the same way, when k = 7r, 
there is no interaction between \n, 0) and the symmetric state so that cu s 
is an eigenvalue of the Hamiltonian which is a zero of the LDOS. By con- 
trast, the coupling with the antisymmetric state leads to two eigenvalues 
-us/ 2 ± ^(u s /2) 2 + Av§ 2 \4) S i\ 2 which yield two poles of the LDOS. 

Consequently, the energy spectrum supports three bands which describe strong 
hybridizations between the three states |n, 0), \a s ) and \a A ). These bands gov- 
ern the dynamics of the lattice which is controlled by a rather large bandwidth 
defined as the energy difference between the high frequency and the low fre- 
quency divergences of the LDOS as Acu = \cos\ + \Ju 2 s + 16t> ^l^sil 2 - The 
numerical analysis reveals that this bandwidth is of about a few times the 
hopping constant <3>. However, as for even v value, the exact diagonalization 
of the full H(k) matrix Eq.(10) reveals that the bandwidth slightly decreases 
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with the number of quanta. It varies from 7.29$ for v — 3 to 6.72$ for v = 51. 
Therefore, this large bandwidth favours a fast decay of the survival probability 
followed by a fast energy transfer. This energy transfer is mediated by bound 
states formed by the superimposition of t>-quanta trapped on the same site and 
v-quanta distributed according to the two states \as) and \oia) which involve 
all the number states \v - 1, 1, 0, 0...), \v - 2, 2, 0, 0...), ...|1, v - 1, 0, 0...). 
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